Simulation of traffic flow at a signalised intersection 
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We have developed a Nagel-Schreckenberg cellular automata model for describing of vehicular 
traffic flow at a single intersection. A set of traffic lights operating either in fixed-time or traffic 
adaptive scheme controls the traffic flow. Closed boundary condition is applied to the streets each of 
which conduct a uni-directional flow. Extensive Monte Carlo simulations are carried out to find the 
model characteristics. In particular, we investigate the dependence of the flows on the signalisation 
parameters. 



I. INTRODUCTION 



Modelling the dynamics of vehicular traffic flow by 
cellular automata has constituted the subject of inten- 
sive research by statistical physics during the past years 
[11) 0, 0, 01 ■ traffic was an earlysimulation target for 
statistical physicists [1, fl 0, 0, i, M, [IH, El, El Q [11 ■ 
Evidently the optimisation of traffic flow at a single in- 
tersection is a preliminary but crucial step to achieve the 
ultimate task of global optimisation in city networks [l^ . 
Nagatani proposed the first model for simulation of two 
crossing streets 17]. Subsequently, Ishibashi and Fukui 
developed similar models [3, [l^. Recently, physicists 
have notably attempted to simulate the traffic flow at in- 
tersections and other traffic designations such as round- 
abouts [13, HH, M, M iM, 25, 26, 27, M M, M M M 
[33I . In principle, the vehicular flow at the intersection of 
two roads can be controlled via two distinctive schemes. 
In the first scheme, the traffic is controlled without traf- 
fic lights [sl,!!!]. In the second scheme, signalised traffic 
lights control the flow. In Ref.[2l|, we have modelled 
the traffic flow at a single intersection with open bound- 
ary conditions on streets. Dependence of waiting times 
on signalisation and inflow rates were investigated in de- 
tails. Along this line of study, our objective in this paper 
is to study in some depth, the characteristics of traffic 
flow and its optimisation in a single intersection with 
closed boundary condition. We will compare our results 
to those obtained in our recent study of a nonsignalised 
intersection [H . In order to capture the basic features of 
this problem, we have constructed a NS cellular automata 
model. This paper has the following layout. In section II, 
the model is introduced and formulated. In sections III 
and IV, the results of Monte Carlo simulations of each 
controlling scheme are exhibited. Concluding remarks 
and discussions ends the paper in section V. 



FIG. 1: Intersection of two uni-directional streets. 



II. 



DESCRIPTION OF THE PROBLEM 



Imagine two perpendicular one dimensional closed 
chains each having L sites. The chains represent ur- 
ban roads accommodating unidirectional vehicular traf- 
flc flows. They intersect each other at the middle sites 
*i = «2 = ■§ on the first and the second chain. With no 
loss of generality we take the flow direction in the flrst 
chain from south to north and in the second chain from 
east to west, (see Fig.l for illustration). Cars are not 
allowed to turn. The discretisation of space is such that 
each car occupies an integer number of cells denoted by 
Lear- The car position is denoted by the location of its 
head cell. Time elapses in discrete steps of At and veloc- 



1,2, 



in which v„ 



ities take discrete values 0, 
is the maximum velocity. 

To be more specific, at each step of time, the sys- 
tem is characterized by the position and velocity con- 
flgurations of cars. The system evolves under the Nagel- 
Schreckenberg (NS) dynamics 
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Let US now specify 
the physical values of our time and space units. The 
length of each car is taken 4.5 metres. Therefore, the 
spatial grid Ax (cell length) equals to m. We take 
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the time step At = 1 s. Furthermore, we adopt a speed- 
Umit of 75 km/h. The value of Vmax is so chosen to give 
the speed-Hmit value 75 km/h or equivalently 21 m/s 
. In this regard, Vmax is given by the integer part of 
21 X Lcar/4.5. For instance, for Lear = 5, Vmax equals 
23. In addition, each discrete increments of velocity sig- 
nifies a value of Au = -^^m/ s which is also equivalent 
to the acceleration. For Lear = 5 this gives the comfort 
acceleration 0.9 Moreover, we take the value of ran- 
dom breaking parameter at p = 0.1. A set of traffic lights 
controls the traffic flow. We discuss two types of signali- 
sation in this paper: fixed-time and traffic responsive. 
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III. FIXED TIME SIGNALISATION OF LIGHTS 

In this scheme the lights periodically turn into red and 
green. The period T, hereafter referred to as cycle time, 
is divided into two phases. In the first phase with du- 
ration Tg, the lights are green for the northward street 
and red for the westward one. In the second phase which 
lasts for T — Tg timsteps the lights change their colour 
i.e.; they become red for the northward and green for the 
westward street. The gap of all cars are update with their 
leader vehicle except those two which are the nearest ap- 
proaching cars to the intersection. These two cars need 
special attention. For these approaching cars gap should 
be adjusted with the signal in its red phase. In this case, 
the gap is defined as the number of cells right after the 
car's head to the intersection point ^. If the head po- 
sition of the approaching car lies in the crossing point 
the gap is zero. We note that at the time step when the 
traffic light goes red for a direction, portion of a car from 
that direction can occupy the crossing point. In this case 
the leading car of the queue in the other direction should 
wait until the passing car from the crossing point com- 
pletely passes the intersection i.e.; its tail cell position 
become larger than . Now all the computational ingre- 
dients for simulation is at our disposal. The streets sizes 
are Li — L2 — 1350 m and we take Lear = 5. The sys- 
tem is update for 2 x 10^ time steps. After transients, two 
streets maintain steady-state currents denoted by Ji and 
J2 which are defined as the number of vehicles passing 
from a fixed location per time step. They are functions of 
global densities pi and p2 and signal times T and Tg. We 
kept P2 fixed in the second street and varied pi. Figure 
(2) exhibits the fundamental diagram of the first street 
i.e.; Ji versus pi. The parameters are specified in the 
caption. 

We observe that for each Tg, Ji linearly increase up 
to pi = p_, then a lengthy plateau region is formed up 
to pi = p+. If Tg increases, p- becomes larger, p+ be- 
comes smaller, and the plateau height increases. This 
seems natural because by increasing Tg the model tends 
to a normal NS model. The emergence of a plateau re- 
gion is associated to defect-like role of the crossing point. 
The asymmetric simple exclusion process (ASEP), as a 
paradigm for non equilibrium processes in one dimension, 
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FIG. 2: Ji vs pi for various values of Tg at fixed p2 = 0.1. 
T = 20 (top), T = 30 (middle) and T = 40 (bottom). 



with one defective site has been investigated within the 
Bethe Ansatz formalism in and matrix product state 
[sj. It is a well-established fact that a local defect can 
affect the low dimensional non-equilibrium systems on a 
global scale. This has been confirmed not only for ASEP 
!34, 38, 39, 40] but also for cellular automata models de- 
scribing vehicular traffic flow [4l|, . After the plateau, 
Ji exhibits linear decrease versus pi in the same manner 
as in the fundamental diagram of a single road. Con- 
cerning the variation of cycle time T, increasing the cycle 
time T gives rise, on an equal basis, to increase both in 



3 



0.25 
^2 0.2 
0.15 
0.1 
0.05 




0.45 
0.4 
0.35 



-T — ^-1 — \ — \ — \ — ^-1 — \ — \ — \ — ^-r~^ — ' — ' — ^~r~^ — ' 

T = 20 _ _□ _ Tg = 5 

^2 V Tg = 15 

I -□- □- B- -□ □- e o e n □ 



A -A^- A- A- A -A-A- A- A -A-i^ A- A 



VVVVVVVVVVVV 



1. I 
V 1 1 



0.2 0.4 0.6 



P2 ~ ^'^ Tg = 20 

Tg = 25 

-o- D-B-H-n □ D-e-o-n-B-a □-□■p 



-A -A- A- & A A A A A--A-A- A A A -A 



VVVVVVVVVVVVVV 



\ \ 
\ I 



0.5 
0.45 

0.4 
0.35 

0.3 

J 0-25 
tot 

0.2 





1 1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 _ 


= -i 


V 


: 

A '' : 


- 


V 


1 - 




□ Tg = 5 

A Tg.10 

V Tg.15 


'\ \ : 
V * 1 

' \ : 


- 


T = 20 


« I 




Pt = 0.1 


\ 










0.2 0.4 0.6 
P, 


0.8 1 



0.4 
0.35 
0.3 

J 0.25 



tot 



0.2 



T = 30 
p,=0.1 



□ Tg = 5 

A Tg.10 

V Tg.15 
--t— - Tg.20 
< Tg = 25 



P, 



FIG. 3: J2 vs pi for various Tg at p2 = 0.1. Top; T=20 and 
bottom: T=30. 



the green and in the red portion of the cycle allocated to 
each street. The results show a notable increase in flows 
when T is increased. This observation does not seem to 
comply to reality. The reason is due to unrealistic nature 
of NS rules. We now consider the flow properties in the 
second street. In Figure (3) we depict the behaviour of 
J2 versus pi for various values of Tg. 

Although p2 remains constant J2 is affected by den- 
sity variations in street 1. For each Tg, J2 as a function 
of pi exhibits two regimes. In the first regime, J2 is al- 
most independent of pi and remains constant up to p+. 
Afterwards in the second regime, J2 exhibits a linearly 
decreasing behaviour towards zero. Analogous to Ji, the 
existence of a wide plateau region indicates that street 2 
can maintain a constant flow capacity for a wide range 
of density variations in the first street. If Tg is increased, 
the green time allocated to the second street decrease so 
we expect J2 to exhibit a diminishing behaviour. This is 
in accordance to what figure (3) shows. In order to find 
a deeper insight, it would be illustrative to look at the 
behaviour of total current Jjot = Ji -I- J2 as a function 
of density in one of the streets. Evidently for optimisa- 
tion of traffic one should maximize the total current Jtot 
therefore it is worth investigating how the total current 
behaves upon variation of the model parameters. Fig. (4) 
sketches this behaviour. 
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FIG. 4: Jtot vs pi for various values of Tg at p2 = 0.1: T = 20 
(Top), r = 30 (middle) and T = 40 (bottom). 



In general, the dependence of total current on Ji de- 
pends on the value of Tg. Except for small values of Tg, 
total current increases with pi then it becomes saturated 
at a lengthy plateau before it starts its linear decrease. 
We have also examined the behaviour of Jtot for other val- 
ues of p2. Figures (5) exhibits the results for p2 = 0.05 
and p2 = 0.5 respectively. 

Our simulations confirm that for small p2 up to 0.1 
total current shows a distinguishable dependence on Tg 
in the entire range of pi especially in intermediate val- 
ues. In contrast, for p2 > 0.1, we observe no significant 
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FIG. 5: Total current Jtot vs pi for various values of Tg at 
T = 30: p = 0.05 (Top) and p = 0.5 (bottom). 



dependence on Tg in the intermediate pi but we observe 
notable dependence for large pi. 



A. Density profile and Queue Formation 

In this section we try to investigate, in brief, the char- 
acteristics of queues which are formed behind the crossing 
point in the course of the red periods. To this end, in the 
following set of figures we exhibit the profile of density 
in the vicinity of intersection. This will provide much 
insight into the problem. In figure (6) we have sketched 
such profiles for three values of green times Tg. 

At each global density pi a high density region, here- 
after referred to as HD, forms behind the crossing point. 
This corresponds to queue formation behind the crossing 
point. The length and height of this HD region grows 
with pi. Keeping pi constant, the length as well as the 
height of the HD region decreases with increasing Tg. 
This is expected since by increasing Tg the cars in the 
first street are given more time to pass the intersection 
and less cars will have to stop at red periods. Right after 
the crossing point, we observe a quick relaxation of the 
density profile to a low density region (hereafter referred 
to as LD). It would be natural to interpret the length 
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FIG. 6: Density profile of the first street for given values of 
parameters which are specified in the figures. Crossing site is 
at 750th cell. 



of the HD region as the average queue length. We note 
that small oscillation in the profiles are not associated to 
inadequacy of simulation time. In fact their emergence is 
related to the updating rules of the NS model. Such os- 
cillations have also been observed in the density-density 
correlation function in the NS model 14311. 
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FIG. 7: Total current vs pi at p2 = 0.05 (top) and p2 = 0.07 
(bottom) for various values of Qc. 



FIG. 8: Jtot vs pi at p2 = 0.15 for various values of Qc- 
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IV. TRAFFIC RESPONSIVE SIGNALISATION 

In this section we present our simulations results for 
the so-called intelligent controllin g sc heme in which the 
traffic light cycle is no longer fixed p^.l45j. In this scheme 
which is some times called traffic responsive, the signali- 
sation of traffic lights is simultaneously adapted to traffic 
status in the vicinity of intersection. This scheme has 
been implemented in simulation of traffic flow at inter- 
sections with open-boundary conditions [13, H^l . There 
exist numerous schemes in which traffic responsive sig- 
nalisation can be prescribed. Three of these schemes are 
discussed in [S^l • Here for brevity we discuss only one of 
these methods. To be precise, we define a cut-off queue 
length Qc- The signal remain red for a street until the 
length of the corresponding queue formed behind the red 
light exceeds the cut-off length Qc- At this moment the 
lights change colour. Apparently due to stochastic na- 
ture of cars movement, the cycle time will be subjected 
to variations and will no longer remain constant. In fig- 
ure (7) we exhibit Jtot versus pi for various values of 
cut-off lengths Qc and p2- 

Analogous to fixed-time scheme, for given p2 a lengthy 
plateau in Jtot forms. The plateau height as well as its 
length show a significant dependence on Qc- higher Qc 



FIG. 9: Variation of cycle lengths in traffic responsive scheme 
for 200 cycles. 



are associated with smaller length and higher current. 
We have also examined larger values of p2- The results 
are qualitatively analogous the above graphs. The no- 
table point is that for p2 larger than 0.1, Jtot do not 
show a significant dependence on p2. Figure (8) depicts 
Jtot versus pi for p2 = 0.15. 

In figure (9) we have depicted the cycle lengths in the 
traffic responsive scheme which vary from cycle to cycle. 
The appearance of such behaviour remarks the adapta- 
tion of traffic lights in the responsive scheme. 

To shed some light onto the problem, we sketch space- 
time plots of vehicles. It is seen that in traffic respon- 
sive scheme, the cars spatial distribution is more homo- 
geneous which is due to randomness in cycle times. 

Lastly, we compare our results to those obtained in 
simulation of a nonsignalised intersection (ssj . In a 
nonsignalised intersection, the cars yield to each other 
upon approaching the crossing point and the priority is 
give to the car which is closer to the crossing point. The 
total current versus pi in a nonsignalised intersection has 
the following form: 

For p2 larger than 0.1, the optimal current is roughly 
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FIG. 10: Space-time plot of vehicles for fixed time (top) and 
traffic responsive schemes (bottom). 
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0.33 in the nonsignalisation scheme whereas in the 
signalised scheme with a fixed time strategy, the optimal 
value of Jtot reaches beyond 0.45 (see Figs. 4 and 5) 
in a fixed time scheme and 0.4 in the traffic responsive 
scheme. It can be concluded that signalisation strategies 
are apparently more efficient in comparison to non- 
signalisation scheme. 

V. SUMMARY AND CONCLUDING REMARKS 

By extensive Monte Carlo simulations, we have investi- 
gated the flow characteristics in a signalised intersection 
via developing a Nagel-Schreckenberg cellular automata 
model. We have considered two types of schemes: fixed- 
time and traffic responsive. In particular, we have ob- 
tained the fundamental diagrams in both streets and the 
dependence of total current on street densities. Our find- 
ings show hindrance of cars upon reaching the red light 
gives rise to formation of plateau regions in the funda- 
mental diagrams. This is reminiscent of the conventional 
role of a single impurity in the one dimensional out of 
equilibrium systems. The existence of wide plateau re- 
gion in the total system current shows the robustness of 
the controlling scheme to the density fluctuations. The 
overall throughput from the intersection shows a signif- 
icant dependence on the cycle time in the fixed time 
scheme and on the queue cut-off length in the responsive 
scheme. Moreover, by plotting the density profiles, we 
have been able to quantify the characteristics of queues. 
Comparison to our previous results for a nonsignalised 
intersection shows the higher efficiency of the signalisa- 
tion strategy. We remark that our approach is open to 
serious challenges. The lack of empirical data prevents us 
to judge how much our results are close to and to what 
degree they differ from reality. Our CA model allows for 
varying space and time grids. By their appropriate ad- 
justing, we are able to reproduce a realistic acceleration. 
We emphasize that our model suffers from non realistic 
behaviour of car movement when the lights turn green. 
The crucial point is to model the braking and accelerat- 
ing as realistic as possible. Empirical data are certainly 
required for this purpose. We expect the system charac- 
teristics undergo substantial changes if realistic yielding 
declaration is taken into account. 



FIG. 11: Jtot vs pi for various values of p2 in a nonsignalised 
intersection. 
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